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Superlattice  (SL)  devices  such  as  infrared  detectors  and  quantum-cascade  lasers  rely  on  efficient  transport 
of  carriers  perpendicular  to  the  SL  layers  by  drift  and/or  diffusion.  While  horizontal  mobilities  are  measured 
routinely,  measurements  of  perpendicular-carrier  mobilities  require  nonstandard  experimental  techniques  such  as 
the  geometric  magneto-resistance.  Here  we  show  how  perpendicular  mobilities  can  be  estimated  from  horizontal 
mobility  measurements  and  calculated  mobilities.  We  treat  low-temperature  horizontal  and  vertical  transport 
in  SL  on  an  equal  footing  by  calculating  both  mobilities  using  the  same  interface  roughness  scattering  (IRS) 
model  from  a  rigorous  solution  of  the  Boltzmann  transport  equation.  The  calculation  is  specialized  to  the  case 
of  InAs/GaSb  SLs,  which  are  of  current  interest  in  the  development  of  third-generation  infrared  detector  focal 
plane  arrays.  The  results  are  compared  to  available  data. 

DOI:  10.1 103/PhysRevB. 84. 155307  PACS  number(s):  73.90.+f,  73.50.Dn,  73.21.Cd 


I.  INTRODUCTION 

The  relationship  between  vertical  and  horizontal  trans¬ 
port  in  superlattices  (SLs)  has  received  renewed  interest 
because  of  type-II  InAs/GaSb  SLs,  which  are  promising 
materials  for  third-generation  infrared  focal-plane  arrays  and 
photodiodes.1-10  The  performance  of  such  devices  requires 
high  vertical  carrier  mobilities  for  efficient  carrier  transport 
through  drift  and/or  diffusion.  However,  vertical  mobilities 
are  not  measured  routinely  since  their  measurement  requires 
nonstandard  and  indirect  experimental  techniques  such  as 
the  geometric  magneto-resistance.11,12  Indeed,  what  is  known 
about  vertical  mobility  is  often  inferred  indirectly  from 
fitting  current- voltage  data.13  Here  we  develop  the  necessary 
formalism  for  the  calculation  of  low-temperature  vertical  and 
horizontal  mobilities  and  specialize  it  to  the  case  of  InAs/GaSb 
SLs.  Then  we  show  how  the  model  results  and  measured 
in-plane  mobilities  can  be  used  to  infer  the  values  of  vertical 
mobilities  for  later  use  in  device  modeling  and  materials 
improvement. 

In  InAs/GaSb  SLs,  hole  transport  in  the  growth  direction 
is  inefficient  since  the  valence  bands  are  relatively  flat  and 
the  corresponding  hole  mass  very  large.14  On  the  other  hand 
the  electron  mass  is  small,  on  the  order  of  the  bulk  mass  in 
InAs  (m  =  0.024mo),  and  electron  motion  is  diffusive,  being 
limited  by  the  usual  scattering  mechanisms  such  as  interface 
roughness  scattering  (IRS),15  phonon  scattering,16  and  point- 
or  extended-defect  scattering.17  Hence,  SL  devices  rely  mainly 
on  electron  transport  for  charge  movement.  Since  vertical  and 
horizontal  electron  transport  in  SLs  is  limited  by  the  same 
scattering  mechanisms,  one  should  calculate  them  using  the 
same  transport  model  so  as  to  compare  them  quantitatively 
and  qualitatively. 

While  there  has  been  a  great  deal  of  work  on  in-plane 
transport  in  quantum  wells,18  there  has  been  no  comparable 
effort  in  SL  transport,  especially  addressing  both  vertical  and 
horizontal  mobilities  and  applying  the  results  to  a  realistic 
system  of  current  interest.19  The  theory  of  diffusive  carrier 
transport  in  SLs  was  developed  by  Mori  and  Ando,20  Dharssi 
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and  Butcher,21  and  others.22,23  Most  theoretical  and  experi¬ 
mental  developments  to  date  concern  horizontal  transport.24 
Here  we  address  ourselves  to  the  issue  of  both  the  vertical  and 
horizontal  transport  in  SLs,  develop  the  required  formalism, 
explore  its  properties,  and  apply  it  to  the  calculation  of  vertical 
and  horizontal  mobilities  in  InAs/GaSb  SLs. 

This  paper  is  organized  as  follows.  In  Sec.  II  the  relevant 
formalism  for  both  the  vertical  and  horizontal  mobilities  in 
SLs  is  derived.  Since  the  relevant  literature  is  scattered  over 
many  publications  or  is  incomplete,  this  paper  will  try  to  be 
self-contained  by  providing  the  salient  detail  in  one  place. 
In  the  case  of  IRS  limited  mobility  it  will  be  shown  that 
in  absence  of  screening,  the  in-plane  relaxation  rates  are 
smaller  than  the  vertical  relaxation  rates;  hence,  in  general, 
if  the  carrier  masses  were  isotropic,  horizontal  mobilities 
would  be  larger  than  vertical  mobilities.  Then,  in  Sec.  Ill 
in-plane  and  vertical-electron  mobilities  are  calculated  for  a 
range  of  SL  and  roughness  parameters  relevant  to  transport 
in  InAs/GaSb  SLs. 


II.  ANALYTICAL  DEVELOPMENT 

Because  of  the  large  hole  masses  of  the  bulk  constituents  of 
the  InAs/GaSb  SL,  the  topmost  heavy  hole  band  in  the  SL  is 
nearly  dispersionless  in  the  growth  direction,  resulting  in  the 
hole  effective  mass  several  orders  of  magnitude  larger  than  the 
electron  mass.14  As  a  result,  holes  are  immobile  and  largely 
confined  to  the  GaSb  layer,  contributing  negligibly  to  the  ver¬ 
tical  current  in  SL  devices  such  as  p-n  or  p-i-n  photodetectors. 
Since  the  photocurrent  in  such  devices  is  carried  mainly  by 
electrons,  the  present  treatment  is  specialized  to  electrons. 

A  transport  calculation  comprises  several  components,  each 
of  which  is  treated  in  turn  subsequently:  the  calculation  of  the 
energy  bands  and  wave  functions — Sec.  A;  the  identification 
and  the  mathematical  treatment  of  the  relevant  scattering 
mechanisms  (here  IRS) — Sec.  B;  the  calculation  of  the  scatter¬ 
ing  rates — Sec.  C;  the  solution  of  the  Boltzmann  equation — 
Sec.  D;  development  of  expressions  for  IRS  relaxation 
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rates — Sec.  E;  development  of  expressions  for  the  mobilities 
and  proof  of  a  theorem  relating  vertical  and  horizontal 
relaxation  rates — Sec.  F;  and  the  calculation  of  relaxation 
rates — Sec.  G. 


A.  Electronic  structure 

The  electronic  structure  of  electrons  in  a  SL  is  found 
from  the  solution  of  the  8x8  EFA  (envelope  function 
Hamiltonian).25  The  SL  consists  of  alternating  layers  of  InAs 
having  width  2 a  and  GaSb  having  width  2b,  with  the  period 
i  =  2a  +  2b.  For  energies  near  the  band  gap  edges,  the  carrier 
wave  function  can  be  approximated  by  a  linear  combination 
of  the  wave  functions  of  the  nearest  bands — the  conduction 
(C),  heavy-hole  (H),  light-hole  (L),  and  spin-orbit  bands 
(SO) — in  each  material,  with  remote  bands  included  via  the 
Lowdin  perturbation  scheme.25  Explicitly,  let  r\\  =  (x,y)  be 
the  in-plane  coordinate  and  z  point  in  the  direction  of  the  SL 
axis.  The  energy  bands  and  wave  functions  are  functions  of 
the  wave  vector  k  =  (k\\,kz),  where  k||  is  the  transverse  wave 
vector,  and  the  perpendicular  wave  vector  —n/t  <  kz  ^  n  / i 
spans  the  Brillouin  zone  in  the  growth  direction.  Overall  the 
envelope  wave  function  for  the  explicitly  included  components 
is  written  as 


<rM,z  |kN,fcz)  =  iAk(rn,z)  = 


1 

Va 


exp(7k,  •  r||) 


8 

x  7 ~2,cn<j)n{kz,z ), 

n=  1 


(1) 


comprising  a  plane  wave  in  the  in-plane  direction  and  a  linear 
combination  of  the  eight-envelope  functions  4>n(kz,z)  in  the 
growth  direction.  Imposition  of  the  usual  boundary  conditions 
gives  rise  to  a  secular  equation  for  energies,  E(k\\,kz),  and 
wave  functions,  {(j)n(kz,z),n  =  1 ...  8},  in  the  8  x  8  envelope- 
function  approximation. 

However,  for  the  system  of  interest  here,  this  level  of  rigor 
is  unnecessary  since  the  calculated  energy  spectra  can  be  well 
approximated  by  the  decoupled  sum26 


h2k2 

E(kn,kz)=  ^A+s(kz),  (2) 

2m\\ 

where  m  n  is  the  effective  mass  in  the  in-plane  direction  and 
s(kz)  is  the  dispersion  in  the  growth  direction.  So,  in  order 
to  make  use  of  this  decoupling  and  to  speed  up  the  transport 
calculation,  the  problem  is  simplified  by  breaking  it  into  two 
parts.  First  the  effective  mass  is  found  by  solving  the  full  8  x 
8  EFA26  at  kz  =  0  for  E(k\\,kz  =  0)  and  fitting  a  parabola  to 
the  calculated  band.  The  parabolic  approximation  holds  very 
well  several  tens  of  meV  above  the  bottom  of  the  conduction 
band. 

The  band  dispersions  e(kz)  in  the  growth  direction  are 
calculated  by  first  setting  kq  =  0  and  neglecting  the  far-off 
bands,  which  decouples  the  8  x  8  EFA  Hamiltonian  into  two 
identical  4x4  Hamiltonians.25  In  turn  each  4x4  Hamiltonian 
splits  into  a  1  x  1  problem  for  the  heavy  holes  and  a  3x3 


Hamiltonian  that  couples  the  C,  LH,  and  SO  bands.27  The  3  x 
3  coupled  EFA  Hamiltonian  reads 


(  Vc(z) 
n/i  Pz 

itlPz 

V3 


VLH(z) 

iAE 

V2 


ijlpz  ^ 
V3 
iAE 

V2 

Vso(z)) 


<  <Pc(z)  \ 
<Plh(z) 
^so(z)  J 


=  £  (kz) 


(  0cfe)  > 
<Plh(z) 
^soiz)  J 


(3) 


where  pz  =  — ih(d/dz )  is  the  momentum  operator;  n  = 
—i(s\px\x) /mo  is  the  Kane  momentum  matrix  element  with 
EP  =  2m0n2;  Vc(z ),  VLH(z ),  and  VSo(z)  denote  the  appro¬ 
priate  strain-dependent,  bulk-band  edges;  and  (pc(z),  <Plh(z ), 
and  ip  so  (z)  are  the  envelope-function  components.  In  terms 
of  perpendicular  and  horizontal  strains,  and  s u ,  and  the 
deformation  potentials  b ,  the  heavy-light  band  splitting  for 
(001)  growth  is  A E  =  2b(s±  —  £||).8,23  In  response  to  strains 
the  Kane  element  changes  to  n  =  n(0)(l  —  s±),  and  the  band 
edges  shift  as  follows.27 


/nx  AQ  1 

VHH  =  vZ+av---AE,  (4) 

A  £2  1 

VLH  =  V^+av—  +  -AE,  (5) 

Vc  =  vP)+ac  — ,  (6) 

Vso  =  V^+av—,  (7) 


where  A£2/£2  =  2s\\  +£±;  ac  and  ay  are  the  hydrostatic 
conduction  and  valence  band  deformation  potentials,  respec¬ 
tively;  A0  is  the  spin-orbit  splitting;  and  V(0)  and  n(0)  denote 
the  unstrained  band  edges  and  the  Kane  matrix  element  in 
the  absence  of  strain,  respectively.  The  InAs/GaSb  band 
overlap — the  difference  between  the  valence-band  edge  of 
GaSb  and  the  conduction-band  edge  of  InAs,  is  taken  to  be 
150  meV.27  The  input  data  used  in  the  model  is  provided  in 
Table  I,  using  data  from  the  model  solid  theory  of  Van  de 
Walle28  and  Landolt-B  or  stein  tables29’  Kane  parameters  n 
are  from  Bastard.25 

Equation  (3)  can  be  simplified  further  to  yield  a 
Schrodinger-like  equation27 


h2  d  Id 
2  dz  m(z)  dz 


+  Vc(z) 


4>c(z)  =  e(kz)</>c(z ), 


(8) 


with  an  “effective”  C/LH/SO  mass  m(z)  defined  by 

,  ,  3[AE2/2-(VLH(z)-e)(VSo(z)-e )] 

m(z)  —  -  \y) 

2n2[2 (Vso(z)  -  A E-e)+  VLH(z)  -  s]  ’ 

with  the  other  two  components  of  the  envelope  wave  function 
given  by, 


ft  (  p  2(VS0(z)  -  s)  -  AE 
0LH  i  yl  3  A£2  -  2 (VLH(z)  -  e)(Vso(z)  -  e) 


x  n 


d^ciz) 

dz 


(10) 
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TABLE  I.  Material  parameters  for  the  model  calculation. 


Ec  (meV) 

Ey  (meY) 

Eso  (meV) 

Yi 

Y2 

Y3 

InAs 

0 

-410 

-790 

-19.67 

-8.37 

-9.29 

GaSb 

960 

157 

-663 

-11.80 

-4.03 

-5.26 

«o(A) 

Cn(10n  dyncm-2) 

C12(10n  dyncm-2) 

C44(10n  dyncm-2) 

InAs 

6.0584 

8.33 

4.53 

3.96 

GaSb 

6.0954 

9.08 

4.13 

4.45 

n2(eV) 

av(eV) 

ac(eV) 

b(eV) 

d(eV) 

InAs 

21.11 

1.0 

-5.08 

-1.80 

-3.60 

GaSb 

22.88 

0.79 

-6.85 

-2.00 

-4.80 

AE  —  (VLH  —  g)  2n  J0c(z) 

*S0(Z)  (AE)2-2(VLH-e)(VS0-e)V3  dz 

(ID 

The  effective  masses  m(z)  [Eq.  (9)]  in  both  materials  are 
shown  in  Fig.  1. 

In  Fig.  1  the  effective  mass  is  strongly  energy  dependent: 
for  example,  for  GaSb  the  mass  becomes  negative  (light-hole 
like)  below  the  valence-band  edge  of  GaSb.  Clearly  it  is 
inaccurate  to  use  bulk  band-edge  masses  in  effective-mass 
like  calculations  far  away  from  band  edges.  Because  of  tensile 
strain  in  In  As,  the  bottom  of  the  conduction  band  in  In  As 
extends  slightly  below  the  bottom  of  the  conduction  band  in 
bulk  In  As. 


The  Schrodinger  equation  for  the  SL  is  solved  subject  to 
the  usual  periodic  boundary  conditions.  In  the  flat-band  case 
the  envelope  function  solutions  are  exponentials  of  the  form25 

<p£{z)  =  CieikAZ  +  C2e~ikAZ  -a^z^a,  (12) 

4>c(z )  =  C?,elkBZ  +  Cj\e~,kBZ  a  ^  z  ^  a  +  2b,  (13) 

where 

kA,B  =  [2 mA,B  (s  -  Vc(z))  /n2]l/2,  (14) 

with  A  =  In  As  and  B  =  GaSb  (Fig.  1).  The  boundary 
conditions  on  the  continuity  of  the  envelope  function  0c  and 
its  derivative  together  with  the  Bloch  condition,  yield 

the  secular  equation 


/  eikAa 

e~ikAa 

_eikBa 

_e~ikBa  \ 

kA  rikAa 

_ kA_p—ikAa 

_  kjLpikBa 

kB  r-ikRa 

c2 

mA 

mA 

mB 

mB 

e~ikAa 

eikAa 

_e~iqd  eikB{a+2b) 

_ e~iqd  e  -i  kB  ( a +2  b) 

c3 

kA  „—ikAa 

i  C 

\mA 

kA  ikAa 
mA  ^ 

_ g—iqd  kB  gikn(a+2b) 

mB 

e~iqd  kB  c-ikn(a+2b ) 
mB  / 

VC4y 

for  the  wave  function  and  the  Kronig-Penney  (KP)  equation 
for  the  energy  eigenvalues,25 

cos kzd  =  cos k&LA  cos  kBEB  —  ^(rj  +  r]~l) 

x  sin&A^A  sin&^L#,  (16) 

where  r\  =  At  the  level  of  approximation  of  the  3x3 
Hamiltonian,  the  LH  and  SO  components  of  the  wave  function 
are  not  continuous.30 

As  an  example,  Fig.  2  shows  several  calculated  band 
structures  for  SLs  with  layer  widths  20A  InAs/40A  GaSb, 
20A  InAs/20A  GaSb,  and  40A  InAs/20A  GaSb  that  span  all 
three  short  (SW)-,  mid  (MW)-,  and  long  (LW) -wavelength 
infrared  (IR)  detection  windows,  respectively.  As  expected, 
the  narrower  the  InAs-layer  width,  the  higher  the  band 
centroid,  and  the  wider  the  period,  the  narrower  the  bandwidth. 
Owing  to  the  band  nonparabolicity  evident  in  Fig.  1,  the 


higher  the  band,  the  heavier  the  carrier  effective  mass;  the 
calculated  numbers  for  the  SWIR,  MWIR,  and  LWIR  SLs  in 
Fig.  2  are  0.0396m0,  0.0297m0,  and  0.0254m0,  respectively. 
Importantly,  a  one-monolayer  variation  of  InAs-layer  width 
changes  the  calculated  bandwidth  only  by  about  7%  for  the 
cases  in  Fig.  1,  which  justifies  our  use  of  the  Boltzmann 
equation  approach.  For  typical  values  of  doping  (1010  cm-2 
per  period),  the  Fermi  levels  are  on  the  order  of  5-30  meV, 
which  translates  into  Fermi  wave  vectors  on  the  order  of  a  few 
tenths  of  the  Brillouin  zone  size.  Therefore,  electrons  occupy 
only  the  ground-state  subband  near  the  bottom  of  the  band. 

Once  the  eigenvalues  of  the  KP  equation  [Eq.  (16)]  are 
located,  the  wave  function  amplitudes  (Ci,C2,C3,C4)  are 
found  from  (15).  The  corresponding  light-  and  spin-orbit 
amplitudes  (Li,L2,L3,L4)  and  (S^S^SA&O  are  calculated 
from  Eqs.  (10)  and  (1 1),  and  then  the  wave  function  normalized 
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as(a)  s  ^  Vlh (GaSb),  kB  =  ikb 

2a(C1C*  +  L\L\  +  SiSl  +  C2C*  +  L2L*  +  S2S*)  +  2  Z?(C3C4*  +  L3L*  +  S3S4*  +  C4C*  +  UL\  +  S4S*)  + 

T  (CiC2  +  L\L/2  H-  *Si ‘S'2  +  C2C*  +  T2T i  sin  2k  ^a / /^4 

+  [(C3c;  +  L3l;  +  S^)^  +  (C4q  +  L4L4  +  S^*)^]  sinh  2  KBb/KB  =  1  (17a) 

and  (b)  e  <  VLH(GaSb) 

2a(CiC;  +  L\L\  +  SiSf  +  C2C*  +  L2L*  +  S2S 2*)  +  2  b(C3c;  +  L3L*  +  S3S 3*  +  C4C4*  +  L4L4  +  5454*) 

+  (CiC|  +  L\L\  +  SiS Z  +  c2c;  +  L2L\  +  S2Si)sin2kAa/kA 

+  [(C3C*  +  L3L*a  +  5354*)easd  +  (C4C*  +  L4l;  +  54^)e-,tB<i]  sin2  kBb/kB  =  1.  (17b) 


In  the  case  of  IRS  the  scattering  rates  depend  on  the  modulus 
of  the  carrier- wave  function  at  the  interfaces  [see  Eqs.  (24), 
(39),  and  (40)  later  in  the  text].  For  the  three  SLs  of  Fig.  2 
the  moduli  of  the  wave  function  at  the  interfaces  are  shown  in 
Fig.  3. 

Figure  3  may  be  understood  as  follows.  The  wave  function 
for  the  ground  state  subband  at  kz  =  0  is  symmetric  with 
respect  to  the  reflection  across  the  midpoint  z  =  0  of  the  InAs 
layer  and  across  the  midpoint  z  =  cl  +  b  of  the  GaSb  layer; 
at  kz  =  7i  /  d  the  wave  function  is  symmetric  about  z  =  0  but 
antisymmetric  about  z  =  a  +  b,  which  is,  therefore,  a  node. 
Hence,  the  value  of  \4>kz(a)\2  must  decrease  as  a  function  of  kz 
from  kz  =  0  to  kz  =  n/d,  as  it  does  in  Fig.  3.  For  the  model 
SLs  the  LH/SO  components  of  the  wave  function  account  for 
less  than  7%  of  the  total  normalization. 

B.  Scattering  mechanism 

Interface  roughness  in  the  form  of  monolayer  fluctuations 
has  been  shown  to  be  the  main  source  of  carrier  scattering 
in  SLs,  affecting  horizontal  transport  in  SLs  up  to  almost 


FIG.  1.  (Color  online)  The  bulk  “effective”  C/LH/SO  mass  m(z) 
Eq.  ((9))  in  the  well  (InAs)  and  barrier  (GaSb)  materials  as  a  function 
of  energy.  The  horizontal  line  is  at  the  value  of  0.024mo  that  is 
appropriate  at  the  conduction  band  edge  of  bulk  InAs.  The  zero  of 
energy  is  at  the  bottom  of  the  conduction  band  of  unstrained  InAs. 
The  strain  conditions  are  those  for  a  20A  InAs/20A  GaSb  SL  grown 
on  a  GaSb  substrate. 


the  room  temperature,  where  polar-optical  phonon  scattering 
becomes  dominant.31  The  characteristic  signature  of  IRS  is  a 
strong  power  dependence  of  mobility  on  layer  width.  Indeed, 
in  the  quantum- well  limit,  Gold’s  theory  predicts  that  mobility 
depends  on  the  sixth  power  of  layer  width.32 

The  fluctuation  of  the  interface  position  at  a  point  rq  is 
characterized  by  height  A(rq)  and  by  the  lateral  correlation 
length  of  interface  fluctuations  A.33  The  fluctuation  gives  rise 
to  a  perturbation  potential  at  each  interface 

SV(t\\,z)  =  Vo  [0(z  -a-  A'(rn))  -  0(z  -  a)] 

+  Vo[6(—z  +  ci  +  A(r  1 1 ))  —  6{— z  +  cl)], 

(18) 

with  the  magnitude  Vo  (given  by  the  height  of  the  InAs 
conduction  well)  and  localized  at  the  well/barrier  z  =  a  and 
barrier/well  interfaces  z  =  —a,  respectively.  To  first  order 
in  the  fluctuation,  the  Taylor  expansion  of  the  perturbation 
yields33 

SV(rn,z)  =  - V0A'(rM )S(z  -  a)  +  V0A(r,| )S(z  +  a), 

(19) 

where  the  fluctuations  A^rq)  and  A(rq)  are  considered 
uncorrelated. 


k/A’1) 


FIG.  2.  (Color  online)  Electron  energy  bands  along  the  growth 
direction  for  three  combinations  of  InAs/GaSb  layer  widths.  For 
typical  carrier  densities,  at  low  temperatures,  the  Fermi  energies  are 
approximately  5-30  meV. 
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In  the  present  paper  the  potential  is  symmetric  (hence 
the  wave-function  moduli  at  z  =  ±a  are  equal),  and  the 
interfaces  are  assumed  to  be  characterized  by  the  same  degree 
of  roughness,  A  =  A',  hence21 


T(k,k)  =  2  x 


2  7T 

h 


V027rA2A2  ,  2  2 

- - - exp(— |k|,  —  k  ||  |2 A2/ 4) 


x  |<Ma)|2|4;(a)|2  M(£(k)  -  £(k'))-  (24) 


Since  |0£/(a)|2  =  |0_^(«)|2,  the  scattering  is 
randomizing,34  T(k\\,kz;kf\\,  —  k'z)  =  T(k\\,kz;k'\\,k'z). 

The  randomizing  property  is  quite  general,  depending  only 
on  the  fact  that  the  autocorrelation  function  of  interface 
roughness  does  not  depend  on  the  vertical  coordinate. 


FIG.  3.  (Color  online)  The  modulus  squared  of  the  electron  wave 
function  at  SL  interfaces  1 0c Ml 2  for  several  combinations  of  layer 
widths  as  a  function  of  the  Brillouin  zone  wave  vector  in  the  growth 
direction  (the  continuous  curve).  Without  the  light-  and  spin-orbit 
wave  function  components,  one  obtains  the  dashed  curves. 


D.  Boltzmann  equation 

When  the  conduction  band  is  wide,  electron  motion  is 
diffusive  and  can  be  described  with  the  Boltzmann  transport 
equation  for  the  electron  distribution  function  /( k).34  In  the 
steady  state  the  Boltzmann  equation  reads35 


The  autocorrelation  function  of  interface  fluctuations  is 
assumed  to  be  Gaussian33  with  the  ensemble  average 

<A(r||)A(r'n»  =  A2exp(-|r,|  -r',||2/A2),  (20) 

and  spectral  density 

S(q)  =  7r  A2  A2  exp(— A2g2/4).  (21) 


C.  Scattering  rates 

The  transition  rate  for  IR-induced  scattering  between 
eigenstates  0k(r||,z)  and  0v(r||,z)  of  the  SL  from  Fermi’s 
Golden  Theorem  is  given  by34 

r(k,k')  =  —  |M(k,k')|2<5(£(k)  -  E( k')),  (22) 

h 

where  the  square  of  transition  matrix  M( k,k')  is 

|M(k,k')|2  =  |(k||,£z|<5V(r||)|k'|,,£'}|2.  (23) 

The  ensemble-averaged  transition  matrix  is  found  to  be 

A  J  exp(/(k'||  —  k,|)  •  (I-,,  —  r'n)) 

x[((A,(r,|)AV||)»l^(fl)|2|^'(a)|2 

+  ((A(r||)A(r'||)))|0iz(— a)|2|0i/(— a)|2]Jr||<ir'||, 

where  ((A(rn)A(r'n)))  is  the  ensemble  average  of  rough¬ 
ness  along  the  z  =  —a  interface  and  ((A^r^A^r'ii)))  for 
the  z  =  a  interface,  the  cross-terms  vanishing  because  of  no 
cross-correlation  between  the  interfaces.  Using  the  Gaussian 
autocorrelation  function,  the  transition  rate  becomes 

n k,k')  =  ^^{A,2A'2exp(-|k||  -k,|||2A'2/4)|^(a)|2 

x  \#kz(a)\2  +  A2A2  exP( — |k||  -  k'n  |2A2/4) 
x  |4,(— a)|2|4z(— a)|2}  S(£(k)  -  £(k')). 


3/  p  1  3/ 

—  =  e¥  • - , 

3 t  c  h  3k 


(25) 


where  electron  acceleration  by  the  applied  electric  field  F  is 
counterbalanced  by  electron  scattering.  The  scattering  term  is 
given  by  the  integral34,35 


=  J [7Xk,k')/(k)(l  -  /( k')) 

-  r(k',k)(l  -  (26) 


where  the  first/second  terms  correspond  to  out-  and  in¬ 
scattering.  In  an  applied  electric  field  F  the  electron- 
distribution  function  can  be  written  as 


f  =  fo  +  g,  (27) 


where  /o  is  the  equilibrium  Fermi-Dirac  distribution  function 
/o  =  [1  +  exp((£  —  EF)/kT)]~l.  In  the  ohmic  regime  the 
departure  from  equilibrium  g  is  expanded  to  first  power  in 
the  applied  field34,35 

e  (  3 E\  3/o 

s  (k)  =  t  (k)  -  •  —  j  A;  (29) 

in  terms  of  the  unknown  relaxation  time  r(k). 

The  relaxation  times  are  derived  by  using  the  microscopic 
reversibility  condition, 


T  (k'.k) 


r  (k  kA  /o(k)(l-/o(kQ) 
1  ’  '  /o(k')(l-/o(k)) 


(30) 


and  inserting  the  distribution  function  [Eq.  (29)]  into  the 
Boltzmann  equation  [Eq.  (26)],  which  yields  the  defining 
equation 


V 


J 


dk'T(k,k’) 


1  ~  /o(k') 
l-/o(k) 


x 


/o(k)(l  -  /o(k)) 

/o(k')(l  -  /o(k')) 


g(k) 

r(k)' 


(31) 
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Hence,  in  general,  the  relaxation  time  is  given  by 


34,35 


1 


T(k) 


(2tt)3 

*r,- 


/* 


T(  k,k) 


1  -  /o(k') 


1  -  /o(k) 
/o(k)(l  -  /0(k))  g(k')' 

/o(k')(l  —  /o(k'))  g(k)J 


(32) 


1.  In-plane  relaxation  rate 
For  the  field  orientation  along  the  x-axis, 
a/0  /i  dE\ 

g(  k)  =  T,,(k  )eFx— - 

7  l|V  7  xdE  \h  dkj 


=  T||(k  )eFx 


/o(l-/o)  /  I  9 E\ 


kT 


nWj’  (33) 


hence  assuming  an  isotropic-parabolic  in-plane  band 
dispersion,34,35 


1 


V 


Tn(k)  (2tt  ): 


X  1 


7 


dk'T(k,kf) 


1  ~  /o(kQ 
1  -  /o(k) 


r||(k')^[i  cos  ip' 


(34) 


T||(k)^||  cos  <p  J 

where  the  angles  between  the  field  direction  and  wave  vectors 
(k||,k'n)  are  (cp,(p'),  respectively. 

Since  cos  cp'  =  cos(p  +  0)  =  cos  <p  cos  0  —  sin  <p  sin  0, 

where  0  is  the  angle  between  vectors  (k||,k'n), 
finally 


1 


V 


T||(k) 


/ 


dk' 


(In)3 

1  -  /o(k') 


,  *ll  Tl|(k )  a 

1  _  1 - ^Tcos0 

k\\  T||(k) 

rc^k'), 


(35) 


1  -  /o(k) 

since  the  sin#  term  drops  out  due  to  parity  considerations. 
Since  IRS  is  elastic,  fo(k')  =  /o(k),  the  relaxation  time 
becomes36 


1 


T||  (kr,kz)  (2^)3 


/ 


kr  {kf^kz) 


where  kr  =  |k||  |. 


x  T(kr,kz;k'r,k';6)k'rdk'rdk'd6,  (36) 


2.  Vertical  relaxation  rate 


For  the  field  applied  in  the  z-direction,  one  similarly 


has 


1 


Tx(k)  (2jt): 


c hf* 


^(kQ  (|f) 
^(k)(||)_ 


Y(k,k);  (37) 


however,  owing  to  randomization,  the  scattering  rate  is  even  in 
k'z  while  3  E/dk'z  is  odd,  so  that  the  in- scattering  term  integrates 
to  zero,  hence34,35 

1  V 


z±  (kr,kz) 


(2  7T)3 

V 

(2^7 


J  T  (k,k')dk'ndk'z 
J  T{kr, kz-,k'r,k'z\6)k'rdk'rdk'zd0. 


E.  Expressions  for  IRS  relaxation  rates 
1.  Vertical  relaxation  rates 

Using  Eq.  (24)  for  the  scattering  rates,  the  integration 
over  the  in-plane  coordinates  k'rdk'rd0'  can  be  performed 
analytically,  yielding23 


1 


r±(E,kz) 


m\\i 


=  2x  A 


f 

Jo 


71 /l 


.2  A  2 


dK\4>K(a)\2h 


A2 


x  exp  (“T  ^  +  k^J  &(E~S  (*0) 

(39) 

where  7o  is  the  zeroth-order,37  modified  Bessel  function  of  the 
first  kind,  kr  =  |k||  |,  and  ©  is  the  Heaveside  step  function. 

2.  Horizontal  relaxation  rates 

Similarly  for  the  horizontal  relaxation  rates  [Eq.  (36)],  one 
must  also  integrate  the  in-scattering  term  containing  cos  Of, 
and  we  derive  the  following  expression, 


1 


r\\(E,kz) 


2x^V02|^(fl)|2j 


n 3 

jt/1 


dk'\<pki(a)\2 


0 


Q2A2 
A  2k’kr 


M^) 


K  qi (^Xk 

kr  T||  (E,kz)  1 


x  exp  (_T  ^  +  k^J  &(E~  8  (*')) 

(40) 

where  I\  is  the  first-order,  modified  Bessel  function  of  the  first 
kind.37  Its  properties  are 

/oW  >  h(x)  >  0, 

with  the  asymptotic  values  for  small  x, 

70  (x)  «  1  +  (x/2)2  H - and  7i  (x)  «  x/2  H - , 

and  for  large  values 


7( 


o  CO  ~  ex /Jinx  ^1  +  ^4 - ^ 


ex /VYnx  I  1 - F  •  ■ 

8x 


and  I\  (x) 


hence  7i(x)  approaches  70(x)  asymptotically  from  below. 

In  the  limit  of  very  large  A,  one  can  use  the  asymptotic 
expressions  and  the  fact  that 


lim 


:  exp 


-(f)2- 


j\  2 


A^oo  2^/jt 

=  8{kr  -  k'r)  oc  8{kz  ±  k'z) 


K) 


to  show  that  r±  1  reaches  an  asymptotic  value  and  r. .  1  grows 


(38)  like  A  2.  In  this  limit  the  scattering  is  purely  specular. 


23 
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F.  Mobilities 

Mobilities  pc  are  found  from  the  calculated  current 

densities,34,35 


j  =  -en/i  •  F  =  —4^3  J  £(k)v(k)dk 

=  -i/r<k)J(F,f)v<k¥k'  (41) 

where  velocity  v(k)  =  ||f  and  carrier  density  n  = 
4^3  f  fo(k)dk.  Hence,  the  two  mobilities  are  found  from 

Mil  =  h!enF\\  and  M±  =  jz/enFz , 


/rx(k)^(ff)Jk 

/  /o(k)Jk 

All  (k)^(afFk 

6  //o(k)Jk 


(42) 

(43) 


The  calculation  will  be  performed  for  temperatures  near 
the  absolute  zero  so  that  dfo/dE  =  —  8(E  —  EF ).  Since  the 
relaxation  rates  themselves  are  temperature  independent, 
an  extension  of  the  calculation  to  higher  temperatures  is 
straightforward,  requiring  only  the  replacement  of  the  delta 
function  8(E  —  EF)  with  dfo/dE  =  /0(  1  —  fo)/kT  and  the 
calculation  of  relaxation  rates  for  energies  within  several  kT  of 
the  chemical  potential  energy;  the  main  challenge  is  numerical. 

Specializing  to  the  case  of  SLs  with  the  energy  band 
dispersion 


h2k2 

E(k)  =  E(kr,kz)=-^+s(kz), 

2m  1 1 

the  carrier  concentration  is 

n  =  _}_  r/l  [Ef  _  e  {k  )]  dk  (44) 

4tt3  h2  J-n/, 

and  the  areal  density  per  period  is  simply  ni .  Similarly,  the 
expressions  for  the  mobilities  are 


1  e 

^X  ~~  n  2tt2^2 


r±(EF,kz) 


1  e  r/l 

Mil  =  ~ ~  ,„2  /  T,|  (EF,kz)[EF  -  s(kz)]dkz. 

n  2jt2h  J-jt/i 


(45) 

(46) 


1.  Theorem  on  mobilities 

From  Eqs.  (39),  and  (40)  the  horizontal  relaxation  rates  are 
seen  to  be  smaller  than  the  vertical  rates,  so  that  the  relaxation 
times  are  related  as 

T||  (E,kr)^T±(EX).  (47) 

If  these  relaxation  rates  were  weighted  equally  in  Eqs.  (42) 
and  (43),  then  in  absence  of  screening  the  IRS  horizontal 
mobilities  would  be  greater  than  the  corresponding  verti¬ 
cal  mobilities  fi\\  ^  ti±.  However,  equal  weighting  requires 
that  the  horizontal  and  vertical  masses  be  equal,  which 


may  not  hold  for  energies  far  from  the  conduction-band 
bottom. 


G.  Calculation  of  the  relaxation  rates 

The  relaxation  rates  are  an  intermediate  result  in  the 
calculation  and,  as  such,  rarely  examined.  However,  in  order 
to  understand  the  calculated  mobilities,  it  is  necessary  to 
examine  the  behavior  of  the  relaxation  rates  as  a  function 
of  SL  and  roughness  parameters.  To  this  end,  the  IRS 
relaxation  rates  were  calculated  for  a  number  of  cases  of 
interest. 


1.  Iteration  experience 

First  we  report  on  the  iteration  process  used  to  solve  the 
integral  equation  for  the  horizontal  relaxation  rates  as  it  bears 
on  the  common  approximations  used  to  calculate  the  in-plane 
mobilities. 

The  iterative  process  is  started  with  the  solution  in  the 
quantum- well  limit,  for  which  k'r  =  kr,  hence  Z\\(E,kr)  = 
Z\\(E,k'r ),  so  the  starting  solution  is 


X  exp  ^  (k'r2  +  k2)^ &(E  -  e'z)  j .  (48) 


The  result  is  then  inserted  in  Eq.  (40)  to  continue  the  itera¬ 
tion  process  by  adding  a  small  fraction  of  the  new  result  to  the 
previous  iteration  until  the  results  are  converged  to  one  part  in  a 
million.  For  A2k2F/2  1,  /i(A2^|/2)  /o(A2&^/2),  so  that 

the  solution  converges  within  a  few  iterations;  for  A2k2F/2 
1 ,  it  can  take  a  hundred  iterations  to  achieve  convergence.  As  an 
example,  Fig.  4  shows  the  iteration  experience  for  the  indicated 
parameters. 

The  zeroth-order  solution  is  seen  to  overestimate  the  final 
result,  especially  near  the  zone  edge. 

Nevertheless  even  with  such  cautious  approach  to  self- 
consistency,  the  next  several  iterations  produce  nonphysical, 
negative  relaxation  rates  in  Fig.  4,  which  points  to  the  failure 
of  the  nonself-consistent  approximation, 


1 


*\\ (E,kz) 


oc 


dk'z\(pk'_(a)\2 


/o 


A2^ 


x  exp 


(49) 


The  final  result  in  Fig.  4  is  physical.  Once  self-consistency 
is  achieved,  however,  further  iterations  can  rapidly  degrade  the 
solution,  leading  to  spurious  oscillations  owing  to  the  round¬ 
off  error. 


2.  Calculated  unscreened  relaxation  rates 

Here  we  examine  the  dependence  of  the  relaxation  rates  on 
SF  parameters.  The  simplest  dependence  is  the  proportionality 
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of  the  relaxation  rates  to  the  square  of  vertical  roughness  A, 
so  that  both  mobilities  are  inversely  proportional  to  A2.  In  all 
cases  examined,  the  vertical  relaxation  rates  are  greater  than 
the  horizontal  rates,  in  agreement  with  the  theorem  derived 
previously. 

Next  we  examine  the  dependence  on  doping  and  the 
correlation  length.  For  typical  values  of  residual  doping,38 
the  Fermi  energy  is  much  less  than  the  bandwidth,  hence  the 
limit  of  the  integration  does  not  extend  to  the  edge  of  the 
Brillouin  zone.  We  distinguish  two  regimes:  A2k2F/2  1  and 

A2k2F/2  1,  cases  (a)  and  (b)  below. 

(a)  Small  values,  A2kF/2  1. 

For  this  case  because  the  in-scattering  term  is  small,  there 
is  a  very  small  difference  between  the  horizontal  and  vertical 
relaxation  rates.  Figure  5  shows  the  results  for  a  20 A  InAs/40A 
GaSb  SL  with  EF  =  25  meV,  A  =  30A,  and  A  =  3 A  so  that 
A2kF/2  =  0.071.  For  small  values  of  x  =  A2kF/2,  I0(x) 

1,  /i(v)  ->  0,  and  exp(— x)  1,  so  that  for  both  mobilities 


1 


1 


r ±(EF,kz)  t||  (EF,kz) 


\ct>kM)\2\K=o(a)\2^.  (50) 


Hence,  the  scattering  rates  are  proportional  to  the 

extent  of  the  available  phase  space  for  scattering  in  the  growth 
direction,  and  the  curves  in  Fig.  5  are  spaced  as  <>JEf.  Since 
the  wave-function  amplitudes  are  largely  independent  of  kz 
for  small  values  of  kz  (Fig.  3),  the  calculated  relaxation  rates 
(Fig.  5)  are  also  relatively  dispersionless. 

(b)  Large  values,  A2kF/2  1. 

For  the  parameters  of  Fig.  6  and  EF  =  25  meV,  A2kF/2  ~ 
5.78.  The  resulting  dependence  on  the  Fermi  energies  and 
vertical  momentum  is  complex  and  best  understood  by 
examining  the  relaxation  rates  near  the  origin  kz=  0  and  at  the 
endpoint  kz  =  kFz,  which  are  plotted  in  Fig.  7.  In  all  cases  even 
for  large  values  of  A,  but  small  Fermi  energies,  the  relaxation 
rates  are  proportional  to  a/^f  [Eq.  (50)]. 

For  larger  Fermi  energies,  one  can  estimate  the  expected 
dependence  as  follows.  Near  the  endpoint  kz  =  kFz  and 


kr  =  0,  so  that  7q(0)  =  1.  Since  the  modulus  |0^(^)l2  Is  a 
relatively  slowly  varying  function,  as  seen  in  Fig.  3,  the  vertical 
relaxation  rate  will  be  proportional  to 

l  r*/i 

—  , — r  oc  e\p(-A2mnEF/2)  /  dk[ 

T±{tF,KFz)  J  o 

x  exp(A2m\\£(k'z)/2)®(EF  —  s(k'z)). 

(51) 

For  the  range  of  doping  considered  here,  the  range  of 
integration  is  limited  by  the  value  of  the  Fermi  wave  vector 
in  the  z-direction,  jt/l  >  kFz  a  \I~E^.  Although  the  integral 
is  nonanalytic,  in  the  limit  of  large  EF  it  is  proportional 
to  exp(A2m\\E F /h2),  so  that  r ±(EF,kFz)~l  ->  constant,  as 
reflected  in  Fig.  7.  In  the  same  limit  the  horizontal  rates  for 
kFz  are  found  to  decrease  faster  than  the  vertical  rates  because 
of  the  growing  importance  of  the  in- scattering  term. 

Near  the  origin  kz  =  0  and  kr  =  kF,  and  for  large  values  of 
Ef ,  using  the  asymptotic  expansion  of  the  Bessel  functions, 
the  vertical  rates  can  be  shown  to  saturate 

T±(EF,kz  =  0)-1  — >  constant, 

and  the  horizontal  rates  to  be  inversely  proportional  to  the 
Fermi  energy 

tii  (efa  =  or1  -*  1  /EP. 

(c)  Dependence  on  layer  widths. 

The  present  theory  reduces  properly  to  the  quantum  well 
limit,  in  which  the  scattering  rates  are  proportional  to  the  sixth 
power  of  layer  width,  here  the  width  of  the  InAs  layer.  The 
analytic  proof  is  given  in  Appendix  A. 

Figures  8  and  9  show  the  calculated  dependence  of  the 
relaxation  rates  on  layer  width  for  two  values  of  correlation 
length.  For  the  short  correlation  length  of  A  =  30  A  (Fig.  8), 
the  vertical  and  horizontal  scattering  rates  are  practically 
indistinguishable  because  the  in-scattering  term  is  negligible. 
For  barrier  widths  2b  =  24,  38,  64  A,  the  rates  can  be 
fitted  to  the  power  laws  r  a  a23  ,a35  ,aAX) ,  respectively,  which 


FIG.  4.  (Color  online)  Iteration  experience  in  solving  the  integral 
equation  for  the  horizontal  relaxation  rates  as  a  function  of  the 
iteration  number  and  the  Brillouin  zone-wave  vector.  The  graph 
highlights  the  inadequacy  of  the  quantum- well  approximation  (zeroth 
iteration)  and  of  the  lack  of  self-consistency  (first  iteration). 


kz  (A-1) 


FIG.  5.  (Color  online)  Calculated,  self-consistent  vertical  and 
horizontal  relaxation  rates  for  a  20A  InAs/40A  GaSb  SL  with  A 
=  30A  and  A  =  3A,  as  a  function  of  the  Brillouin  zone  wave  vector 
and  Fermi  energy.  The  curves  terminate  at  the  value  of  kz,  for  which 
Ef  =  s(kz). 
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FIG.  6.  (Color  online)  Calculated  vertical  and  (self-consistent) 
horizontal  relaxation  rates  for  a  20A  InAs/40A  GaSb  SL  with  A  = 
200A  and  A  =  3A,  as  a  function  of  the  Brillouin  zone  wave  vector 
and  Fermi  energy. 


approaches  the  quantum- well  limit  as  InAs  wells  become  more 
separated.  Using  only  the  data  for  InAs  widths  greater  than 
30  A,  the  fitted  power  laws  r  a  a2  6, a42, a5  0  approach  the 
quantum- well  limit  even  closer. 

For  longer  correlations  lengths  A  =  200  A  (Fig.  9),  the  ver¬ 
tical  (t_l  oc  a23  ,a33  ,a40)  and  horizontal  (t|  |  a  a2  3 , a3  6 , a3  9) 
scattering  rates  differ  markedly,  but  their  power  dependences 
are  similar  to  each  other  and  to  those  for  the  shorter  correlation 
length  (Fig.  8).  Using  only  the  data  for  InAs  widths  greater  than 
30  A,  the  fits  are  r±  oc  a26, a42, a5  0  and  t\\  oc  a26, a43, a51. 
The  similarity  of  these  power  laws  underscores  that  the  InAs- 
width  dependence  of  the  scattering  rates  derives  primarily  from 
the  wave-function  amplitudes. 


FIG.  7.  (Color  online)  Calculated  vertical  and  horizontal  relax¬ 
ation  rates  for  a  20A  InAs/40A  GaSb  SL  with  A  =  200A  and  A  = 
3A,  as  a  function  of  Fermi  energy  near  the  center  of  the  Brillouin 
zone  and  at  the  endpoint  kFz. 
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FIG.  8.  (Color  online)  The  calculated  electron  relaxation  rates  for 
several  GaSb  barrier  widths  as  a  function  of  InAs  well  width  for  a 
short  correlation  length:  A  =  30  A,  EF  =  25  meV,  A  =  3 A,  and 
kz  =  0. 

III.  CALCULATED  MOBILITIES 
A.  Unscreened  results 

A  number  of  mobility  curves  were  calculated  as  a  function 
of  SL  parameters,  using  converged  relaxation  rates.  For  this 
section  only  the  parallel  mass  was  fixed  at  m\\  =  0.024mo. 
Figure  10  shows  the  calculated  mobilities  as  a  function  of 
InAs-layer  width  for  correlation  lengths  A  =  30  and  200  A  and 
Ef  =  3  meV.  Since  with  increasing  layer  width  the  interface 
is  an  increasingly  smaller  fraction  of  the  total  width,  the 
calculated  mobilities  increase  with  InAs-layer  width,  which  is 
the  hallmark  of  IRS.32  For  InAs  widths  greater  than  60  A,  the 
conduction-band  level  falls  below  the  top  of  the  bulk- valence 
band  in  GaSb,  so  that  the  electron  wave  function  no  longer 
decays  exponentially  but  oscillates  in  GaSb.  In  the  3  x  3  EFA 
model  the  transition  from  the  decaying  to  oscillatory  behavior 
inside  the  GaSb  layer  is  smooth,39  with  mobilities  rising  at 


InAs  (A) 


FIG.  9.  (Color  online)  The  calculated  electron  relaxation  rates  for 
several  GaSb  barrier  widths  as  a  function  of  InAs  well  width  for  a 
long  correlation  length:  A  =  200  A,  EF  =  25  meV,  A  =  3 A,  and 
kz  =  0. 
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FIG.  10.  (Color  online)  Calculated  vertical  and  horizontal  un¬ 
screened  mobilities  as  a  function  of  InAs-layer  width  for  the  indicated 
parameters;  the  parallel  mass  m\\  =  0.024m0. 


an  ever  increasing  rate  as  the  In  As  layer  widens.  For  A  = 
30  A  the  vertical  mobilities  become  slightly  larger  than  the 
horizontal  mobilities  because  of  the  unequal  weighting  of  the 
corresponding  relaxation  rates  by  the  velocities  in  the  mobility 
integral. 

On  the  other  hand  the  calculated  mobilities  decrease  with 
GaSb  layer  width  (Fig.  11)  since  the  electron  wave  function 
cannot  penetrate  into  the  barrier  material  and  so  is  more 
susceptible  to  InAs-layer  width  fluctuations. 


B.  Screened  results 

In  this  section  we  use  the  data  for  the  four  samples  we 
studied  earlier.  The  experimental  sample  parameters  such  as 
layer  widths  and  horizontal  mobilities  are  listed  in  Table  II.38 
In  addition  we  calculate  the  Fermi  energies  using  the  bands 
calculated  with  the  3x3  model,  and  the  parallel  masses  and 
band  gaps  are  calculated  using  the  full  8x8  EFA  model.  Since 
we  are  dealing  with  n- type  transport,  we  do  not  include  the 
p-type  sample,  labeled  as  SL1  in  Ref.  38,  whose  narrow 
width  and  measured  low  mobility  led  to  the  high-power 
dependence  of  mobility  on  layer  width  reported  in  Ref.  38.  Its 
inclusion  at  the  time  was  motivated  by  the  fact  that  we  were 
using  Gold’s  theory  that  is  appropriate  to  isolated  quantum 
wells,  which  moreover  does  not  distinguish  between  n-  and 
/7-type  conduction  and  is  independent  of  parallel  mass.  The 
challenging  /7-type  transport  merits  an  independent  study. 


FIG.  11.  (Color  online)  Calculated  vertical  and  horizontal  un¬ 
screened  mobilities  as  a  function  of  GaSb-layer  width  for  the  indicated 
parameters;  parallel  mass  m\\  =  0.024 m0. 


For  comparison  with  data  it  is  necessary  to  take  into  account 
the  electron  screening  of  IRS.  Here,  we  adapt  the  Thomas- 
Fermi  form  of  screening,  which  requires  that  the  transition 
matrix  be  multiplied  by  the  screening  factor40 

\q  +  qs )  \|k||  —  k  |||  +  qs  J 

=  /  VWjz  2krKC0se'  +  W  V  (52) 

\  \[W~  2krK  C0S  0'  +  K2  +  9s) 

where  qs  is  the  Thomas-Fermi  screening  wave  vector  qs  = 

2 

and  s  is  the  bulk  dielectric  constant.  For  the  numbers 
appropriate  to  In  As,  the  Thomas-Fermi  screening  length  is 
about  q~x  ~  180  A.  Given  the  carrier  densities  in  Table  II, 
( kp)~l  ~  200A,  which  is  on  the  order  of  q~l ,  so  that  electronic 
screening  cannot  be  neglected.  Inclusion  of  screening  makes 
the  angular  integrals  for  the  scattering  rates  nonanalytic. 

The  overall  effect  of  screening  is  to  reduce  the  scat¬ 
tering  rates  and  enhance  mobilities.  However,  the  verti¬ 
cal  and  horizontal  scattering  rates  are  affected  differently 
because  of  the  angular  dependence  of  Sc(q ).  In  fact  the 
vertical  rates  are  reduced  more  since  screening  is  largest 
for  backward  scattering,  0r  =  n,  where  1  —  (k'r / kr)  cos  0'  = 
1  +  (k'r/kr)  so  that  horizontal  rates  are  accentuated  with 


TABLE  II.  Data  for  the  four  samples  in  this  study,  including  the  experimental  layer  widths  and  horizontal  low-temperature  (10  K)  mobilities, 
as  well  as  calculated  Fermi  energies,  parallel  masses,  and  band  gaps  from  the  8  x  8  EFA  model. 


Superlattice 

InAs 

width  (A) 

GaSb 
width  (A) 

Sheet  carrier 
concentration  per  period 

(xlO10  cm-2) 

Fermi 

energy  (meV) 

Calculated 

parallel 

mass  (8x8  EFA) 

Measured 

mobilities 
(cm2 /V- sec) 

Calculated 
band  gaps 
(meV)  InSb-IFs 

SL2 

41.1 

24.0 

1.8 

13.38 

0.0306 

747 

120.2 

SL3 

53.6 

24.0 

2.4 

15.19 

0.0291 

5043 

66.6 

SL4 

62.3 

24.0 

3.2 

17.73 

0.0279 

8980 

47.3 

SL5 

72.7 

24.0 

3.0 

16.62 

0.0261 

15450 

43.0 
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FIG.  12.  (Color  online)  The  calculated  vertical  (continuous)  and 
horizontal  (dashed)  mobilities  for  four  samples  for  which  low- 
temperature  in-plane  mobilities  (denoted  by  squares)  are  available 
(see  Table  II). 

respect  to  vertical  rates.  Hence,  with  screening  included,  it 
is  possible  for  vertical  mobilities  to  be  larger  than  horizontal 
mobilities. 

As  an  example  of  a  screened  calculation,  Fig.  12  shows 
the  calculated  mobilites  for  the  parameters  of  the  four  samples 
studied  earlier  (Table  II).38  Given  the  experimentally  measured 
carrier  concentrations,  the  Fermi  levels  for  all  the  samples  were 
calculated  from  Eq.  (24)  using  the  theoretically  calculated 
band  dispersions  s(kz).  For  the  sheet  carrier  concentrations 
per  period  on  the  order  of  3.0  x  1010  cm  2  (Table  II),  the 
Fermi  energies  for  the  four  samples  are  between  13.38  and 
17.73  meV.  Vertical  roughness  A  is  taken  to  be  one  monolayer, 
and  the  correlation  length  is  left  as  a  free  parameter.  Figure  12 
spans  the  range  of  correlation  lengths  found  experimentally  by 
scanning  tunneling  microscopy  by  Lew  et  al.  41  (A  =  1.9  ± 
0.1  A,  A  =  1 12  ±  16  A;  A  =  2.8  ±  0.2  A,  A  =  174  ±  21  A)  for 
InAs  on  InGaSb  and  (A  =  3.2  zb  0.2  A,  A  =  301  d=  39  A;  A  = 
4.3  ±  0.2  A,  A  =  327  d=  38  A)  for  InGaSb  on  InAs  interfaces 
in  (1 10)/(1— 10)  cross-sections,  and  by  Feenstra  et  al.  42  (A  = 
1.8  A,  A  =  19  A)  for  InAs  on  GaSb  and  (A  =  2.5  A,  A  = 
40  A)  for  GaSb  on  InAs  interfaces.  The  calculated  mobilities 
together  with  the  measured  data  are  shown  in  Fig.  12. 

As  in  the  case  of  quantum  wells,32  the  horizontal  mobility  is 
a  double- valued  function  of  A,  with  a  minimum  corresponding 
to  the  value  of  A  comparable  to  the  Fermi  wavelength,  where 
the  carriers  experience  the  greatest  scattering.  To  the  left 
of  the  minimum,  the  electron  averages  over  the  small-scale 
roughness,  so  the  mobility  is  high.  To  the  right  of  the  minimum, 
the  electron  wavelength  is  smaller  than  A,  and  the  interface 
appears  smooth  for  in-plane  transport,  and  the  mobility  is 
high  again.  Vertical  mobility  reaches  a  steady-state  value, 
since  once  the  interface  is  sufficiently  smooth,  scattering 
becomes  specular  and  depends  only  on  vertical  roughness.23 
For  low  values  of  A,  the  vertical  mobilities  are  larger  than  the 
horizontal  mobilities  by  a  factor  of  approximately  3  to  2.  For 
A  around  225  A,  the  horizontal  mobilities  are  twice  as  large  as 
the  vertical  mobilities;  for  A  around  300  A,  that  ratio  is  about 
3  to  1. 


FIG.  13.  (Color  online)  Calculated  IRS-limited  mobilities  for  a 
41.1  A  In  As/24  A  GaSb  SL  with  screening  for  two  values  of  the 
autocorrelation  length  as  a  function  of  the  Fermi  level.  The  carrier 
concentrations  corresponding  to  the  given  Fermi  level  can  be  read  off 
along  the  right-hand  scale. 


Using  the  assumed  value  of  A,  two  values  of  A  will 
fit  the  experimental  mobilities  for  the  three  widest  samples, 
but  no  value  of  A  will  fit  the  experimental  mobility  for  the 
narrowest  SL.  Since  mobilities  are  inversely  proportional  to 
the  square  of  vertical  roughness,  one  can  achieve  agreement 
with  experiment  by  simply  adjusting  the  value  of  A.  However, 
we  also  considered  possible  problems  with  the  measurement 
and  theory.  Experimentally,  we  find  that  at  InAs  widths  greater 
than  40  A  samples  type-convert  from  p  to  n ,  which  may 
indicate  a  more  complicated  physics  than  assumed  in  the 
theory.  In  addition,  there  is  great  sample-to- sample  variability 
in  mobilities,  depending  on  the  exact  growth  conditions. 
Theoretically,  there  may  also  be  other  sources  of  scattering.  For 
example,  we  considered  the  effect  of  possible  alloy  scattering 
due  to  antimony  interdiffusion  into  InAs43  but  found  the 
strength  of  such  scattering  to  be  two  orders  of  magnitude 
too  small  to  account  for  the  observed  discrepancy  between 
experiment  and  theory  (see  Appendix  B). 

Since  the  ratio  of  the  two  mobilities  is  independent  of  A, 
curves  like  those  in  Fig.  12  make  it  possible  to  estimate  the 
values  of  vertical  mobilities  given  the  measured  value  of  the 
corresponding  horizontal  mobility,  if  an  additional  piece  of 
information  such  as  horizontal  roughness  or  correlation  length 
are  known.  Alternatively,  horizontal  mobility  measurements 
over  an  extended  temperature  range  may  provide  enough 
additional  information  to  make  an  independent  determination 
of  A  unnecessary. 

An  example  of  the  concentration  dependence  of  mobilities 
is  given  in  Fig.  13  for  a  4 1.1  A  InAs/24A  GaSb  SL  for  two 
values  of  correlation  length.  For  a  small  correlation  length,  A  = 
30  A,  both  mobilities  decrease  as  a  function  of  the  Fermi  level 
or,  equivalently,  the  carrier  density.  For  the  long  correlation 
length  of  A  =  200  A,  the  horizontal  mobility  first  decreases 
and  then  rises  again,  while  the  vertical  mobility  decreases  to 
an  asymptotic  value.  This  behavior  is  consistent  with  that  in 
Fig.  12,  except  that  in  Fig.  13  the  correlation  length  is  fixed  and 
the  minimum  is  found  at  the  Fermi  energy  corresponding  to  the 
Fermi  wavelength  on  the  order  of  the  fixed  correlation  length. 
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For  A  =  30A  the  minimum  takes  place  at  carrier  densities  Inserting  this  dependence  into  the  expression  for  the 

beyond  the  range  of  the  graph.  relaxation  rate  gives 


IV.  CONCLUSIONS 

The  theory  of  low-temperature  electron  transport  in  SLs 
was  presented,  treating  vertical  and  horizontal  mobilities  on 
an  equal  footing  within  the  same  IRS  model.  Low  temperature 
vertical  and  horizontal  electron  mobilities  in  InAs/GaSb  SLs 
were  calculated  by  solving  the  corresponding  Boltzmann 
equations  using  the  energy  spectra  and  wave  functions  ob¬ 
tained  from  the  solution  of  the  coupled  envelope-function 
Hamiltonians  (3x3  for  the  vertical  and  8  x  8  for  the  horizontal 
dispersions)  for  the  SL.  The  behavior  of  the  scattering  rates 
and  mobilities  was  examined  as  a  function  of  SL  parameters 
and  explained  in  terms  of  the  underlying  physics.  In  absence 
of  screening  a  theorem  was  derived  showing  that  vertical 
scattering  rates  are  greater  than  the  horizontal  rates  so  that  the 
vertical  mobilities  should  be  generally  smaller  than  horizontal 
mobilities.  With  screening  included,  backward  scattering  in 
horizontal  transport  is  accentuated,  so  that  no  general  predic¬ 
tions  can  be  made.  The  calculated  mobilities  were  found  to  be 
strong  functions  of  the  interface  roughness  parameters,  carrier 
concentration,  and  carrier  screening.  Theoretically  calculated 
mobilities  compared  favorably  with  limited  low-temperature, 
horizontal  mobility  data  on  four  selected  samples.  Using  the 
present  results  and  the  measured  values  of  horizontal  mobil¬ 
ities,  it  is  possible  to  estimate  the  value  of  the  corresponding 
vertical  mobilities. 
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APPENDIX  A 


_1_  =  (-_<_» WAn  [■" 

T||(k)  \Anh{m*)2{2af  )  J-n/1 

x  J  (1  —  cos0')exp(— A2(kj|  —  k||)2/4) 
x  8(s  —  e'y2kj|,  (A4) 

and  using  e(k||)  =  h2k^/2m*,  the  k-  integral  gives  2n/l. 
Integrating  over kj|  withi5(s  —  s')  =  m*S(k^  —  k\\)  /  (fi2k\\),  the 
final  result  is 

1  / TC4hA2A2\  f  .  A 

T\\(k\\)  \2m*(2af)  J  V  7 

(o  9 1  —  COS  0'\ 

— A2k2 - - - j  dO  ,  (A5) 

so  that  T\\  and  mobility  are  proportional  to  the  sixth  power 
of  layer  width.  This  agrees  with  the  expression  for  quantum 
wells  derived  by  Gold.32  The  relaxation  rate  in  the  vertical 
direction  has  the  same  functional  dependence  but  without  the 
cos  O'  term.  Hence,  the  present  theory  correctly  reduces  to  the 
quantum- well  limit  as  the  barrier  height  becomes  infinite. 

APPENDIX  B:  ESTIMATE  OF  ALLOY  SCATTERING 

In  this  section  we  develop  an  expression  for  the  relaxation 
time  in  SLs  for  alloy  scattering  in  order  to  compare  the 
strengths  of  alloy  and  IRS  scatterings.  In  the  binary-binary  SL, 
alloy  scattering  can  be  due  to  Sb  interdiffusion  into  In  As  layers, 
which  results  in  an  InAsi_xSbx  alloy  with  the  alloy  profile  x(z). 
In  order  to  calculate  the  relaxation  rates  for  alloy  scattering,44 
the  potential  in  the  current-carrying  region  is  written  in  the 
virtual  crystal  form  as45 


In  order  to  evaluate  the  quantum- well  limit  of  the  expression 
for  mobility,  it  is  necessary  to  find  the  product  of  well  height 
and  the  fourth  power  of  the  modulus  of  the  wave  function  in 
the  limit  of  infinite  barrier  height.21  Let  the  quantum  well  of 
height  Vo  be  located  between  —a^z^a.  The  electron-wave 
function  in  the  ground  state  is 

I  A  cos  kz  —a^z^a 

Be~KZ  |z|  >  a  ’  (A1) 


Vo(r)  =  J2[X  VlnSb(r  -  a)  +  (!  -  xWlnAs(r  -  a)],  (Bl) 

a 

where  the  sum  is  over  all  the  molecular  sites.  Scattering  is  due 
to  random  deviations  from  the  average  crystal, 

5  Vb(r)  =  X  ca[ymsb(r  -  a)  -  VInAs(r  -  a)],  (B2) 

a 


where  k  =  y/2m*E  /h2  and  k  =  y/ 2m*(Vo  —  E)/Ti2.  The 
boundary  conditions  on  the  wave  function  and  its  derivative 
furnish  the  relation  (p(a)/cp'(a)  =  — 1/k ,  so  that 


VoV(«)  =  V2  (<p'  (a)/K)4  =  V2 


E) 


x  ( kA  sin  kay . 


(A2) 


In  the  limit  of  infinite  barrier  height,  the  amplitude  A  = 
1  /y/a,  k  =  n/2a,  and  sin  £<2  =  1,  so  that 


lim  V^(pA{a)  =  (h27t2 / m*)2 /(2a)6 .  (A3) 

Vo— ^00 


where  ca  is  a  random  function  equal  to 


C;i  — 


l  —  x  for  an  InSb  molecule  at  a, 
— x  for  an  InAs  molecule  at  a, 


(B3) 


and  the  scattering  potential  is 


VinSb(r  -  a)  -  V inAs(r  -  a)  =  v0w05(r  -  a),  (B4) 

where  coo  =  a3/ 4  is  the  volume  per  anion-cation  pair46  and  Vo 
is  the  conduction-band  offset  between  InSb  and  InAs. 
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The  matrix  element  of  the  scattering  potential  is  calculated 
to  be 


|M(k,k')|2  = 


(^r)  /  [1-^)]*(z)i<M*)i2 


X  \(pk'Xz)\2dz, 


(B5) 


which,  like  IRS,  is  also  randomizing  and  elastic.  Since  the 
scattering  mechanism  is  isotropic  and  independent  of  the 
scattering  angle  O',  the  in- scattering  term  for  the  horizontal 
relaxation  rate  integrates  to  zero,  so  that  the  horizontal  and 
vertical  relaxation  times  are  equal, 

1  1  im\\ 

Tj_(E,kz)  T||  (E,kz)  2nh3 

x  J  l(kz,k'z)@(E  -  e'z)dk'v  (B6) 

where 


(1  -  x{z))x(z)\(t>kz{z)\2\<t>k[{z)\2dz 


(B7) 


is  the  overlap  of  the  grading  profile  with  the  initial  and  final 
state  wave  functions.  Equation  (B6)  agrees  with  the  similar 
expressions  derived  earlier  for  quantum  wells.47-49 


To  compare  alloy  to  IRS  scattering,  we  assume  that  Sb 
diffuses  into  InAs  at  both  interfaces  by  about  one  monolayer, 
8  ~  3A,  so  that 


x(z)  = 


1  -(z  +  a)/8 
0 

1  +  (z-a)/8 


—a^z^—a  +  8 
—a  +  8^z^a  —  8' 
a  —  8  ^  z  ^  a 


(B8) 


Inserting  this  profile  into  Eq.  (B6),  one  may  approximate 
the  relaxation  rate  by 


1 

t±(E,kz) 


/w||£ 

V  n3 


vl\<pk,(a)\ 


8co  o 
3i x 


x  /  |<fe(a)|2©(£  -  e'z)dk’z. 

Jo 

(B9) 

The  difference  between  Eq.  (B4)  for  alloy  scattering  and 
Eqs.  (39)  and  (40)  for  IRS  scattering  is  mainly  in  the  strength 
parameters:  Vq8coq/37T  for  alloy  scattering  and  2VqQ2A2  for 
IRS.  Inserting  typical  numbers  for  the  InAs/GaSb  SL  shows 
that  alloy  scattering  is  about  100  times  weaker  than  IRS 
scattering. 
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